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Abstract 

The optimized (5-expansion is used to study vacuum polarization effects in the 
Walecka model. The optimized 5-expansion is a nonperturbative approach for 
field theoretic models which combines the techniques of perturbation theory and 
the variational principle. Vacuum effects on self-energies and the energy density of 
nuclear matter are studied up to 0{5'^). When exchange diagrams are neglected, 
the traditional Relativistic Hartree Approximation (RHA) results are exactly re- 
produced and, using the same set of parameters that saturate nuclear matter in 
the RHA, a new stable, tightly bound state at high density is found. 
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Although hadrons are not elementary particles, the use of renormalizable relativistic 
quantum field models employing pointlike hadrons, such as the Walecka model [|l|, for 
studying the properties of hadronic matter at high density and/or temperature is im- 
portant in several aspects. Perhaps one of the most important aspects is the necessity 
of understanding the vacuum of such models, since one expects that the description of 
highly excited matter in terms of such models must breakdown at some scale and then 
quark-gluon degrees of freedom must be invoked for a proper treatment of the system. 
Quantum fiuctuations in the Walecka model have been studied at the Hartree level 0, 
but severe difficulties arise when nonperturbative exchange diagrams (Fock diagrams) 
are considered. Although Bielajew and Serot |^ have set up an appropriate framework 
for the renormalization of the Hartree-Fock (HF) equations, no explicit calculation em- 
ploying this has been performed so far. Bearing in mind the difficulties associated with 
the renormalization within the HF approximation we give in this Letter the first step to- 
wards the implementation of a different nonperturbative approach for studying vacuum 
fluctuations in hadronic models. Namely, we use the 6 expansion [^]; more speciflcally 
we use the optimized linear 6- expansion 0. 

The standard application of the linear 5-expansion p, |^, ^ to a theory with action 
S starts with an interpolation deflned by 

S{S) = (1 - 6)So{fi) + SS = Soifi) + S[S - Soifi)], (1) 

where S'o(/i) is the action of a solvable theory. The action S{6) interpolates between the 
solvable So{fi) (when 6 = 0) and the original S (when 6 = 1). Since Sq is quadratic in 
the flelds, arbitrary parameters (/i) with mass dimensions are required for dimensional 
balance. The evaluation of a physical quantity P is performed by considering the term 
6[S — 5'o(yu)] as a perturbation whose order is labeled by 6 which is set to unity at the 
end. In practice, the perturbative expansion in powers of 6 will be truncated at a given 
order implying that the quantity P will have a residual dependence on the unknown 
parameters. Since /i does not belong to the original theory it will have to be flxed 
according to some criterion and different methods have been proposed. Among them 
the Principle of Minimal Sensitivity (PMS) ^ offers a particularly attractive way of 
optimizing the theory by requiring P(/i) to be evaluated at the point where it is less 

2 



sensitive to small variations of /i: dP{fi)/dfi = 0. With this variational procedure /i 
becomes a function of the original parameters of the theory yielding non-perturbative 
results. 

The different forms of the 6 expansion have been successfully applied to many dif- 
ferent problems in quantum mechanics, particle theory, statistical physics and lattice 
field theory, and its convergence has been recently proved for quantum mechanical prob- 
lems More recently, the optimized linear 5-expansion was used for applications in 
$1, QED and Yang-Mills theories [|10|. In a recent paper ^ it was demonstrated that in 
truncating the 6 expansion at and neglecting vacuum effects, one can readily re- 

produce the results of the standard self-consistent Dirac-Hartree-Fock approximation 
(in this approximation vacuum effects are also neglected) for the equation of state of nu- 
clear matter. In this Letter we go one step further by including vacuum effects to (9(5^) 
and demonstrate that one can reproduce the standard relativistic Hartree results (RHA). 
In addition, from the PMS applied to the renormalized energy density of nuclear matter, 
a new stable state is found at high density. One of the stengths of the 5-expansion is 
that one can proceed to beyond leading order with considerably less computational effort 
than in the traditional HF approximation since in one deals in this approach with a finite 
number of Feynman graphs as implied by the perturbative nature of the calculation. 

We start with the Lagrangian density of the Walecka model [jl], ||: 

£ = ^{t ^-M + g,<f)-g,YV^)^ + ^{d^(l)d^<f)-ml<j)^) 

-\f^uF^'' + ]^mlV,V^ + V) + £cT, (2) 

where if) represents the nucleon field operator, and are respectively the field oper- 
ators of the scalar and vector meson and F^^ = d^V^ — dyV^. The massive vector field 
is coupled to a conserved baryon current, rendering the model renormalizable in 3+1 
dimensions. In order to minimize many body effects the term V), which describes 
mesonic self interactions, was set to zero in the original work of Walecka and we shall 
keep the same convention here. The Lagrangian density Cct contains all the countert- 
erms needed to render the model finite which, for the purposes of the present paper, is 
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given by: 



n=l 



(3) 



We are primarily interested in studying effects of vacuum fluctuations in nuclear matter. 
In particular, we will consider the energy per nucleon, which is related to the energy 
density by: 

S = l-j£x (<^|T°°|^> - <vac|T°°|vac>) =£^ + + (4) 

where |\E'> is the interacting ground-state of nuclear matter, |vac> is the vacuum state 
(zero density), includes the piece corresponding to Cct- The terms S^, S'^ and 
refer respectively to baryon, scalar-meson, and vector-meson contributions, which are 
given by: 



{Tr [7°fc° 
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where Am(A;),m = s,v are the Fourier transforms of Green's functions of the Klein- 
Gordon operator. It is important to note that we are not using the nucleon equation of 
motion and, therefore, the above expressions differ from the usual ones . The reasons 
for not using the nucleon equation of motion will be discussed below. The terms £ct'^ ^ill 
be considered in renormalization of the energy density. In obtaining these expressions, 
we have eliminated the meson field operators and in favor of the nucleon field 
operators by integrating the meson Euler-Lagrange equations. 

The problem of calculating the energy density consists therefore in finding the in- 
teracting nucleon propagator in medium. The first step towards the calculation of the 
propagator with the 5 expansion is to define the interpolated Lagrangian. According to 
Eq. (p, one has: 

C{5) = {l-5)C, + 5C = C,+ 6tP{g,<p - g,YV, + /i)^/' + Cct{S). (6) 



where we have chosen Cq to be: 

where f2 = M + /i. Notice that the (5-expansion interpolation could also have been done 
in the mesonic sector. However, we have chosen to eliminate the mesonic fields by in- 
tegrating their equations of motion; in this way mesonic self-energies are automatically 
taken into account. Of course, the integration in favor of the nucleon fields can be done 
only in the absence of mesonic self-interactions [Note that since we have set V) =Q 
the mesonic self-interactions appear in Cct only]. This leaves us with only one un- 
known parameter, /z, which will be fixed by the PMS condition applied to the energy 
density. In short, one is performing a variational calculation of the energy density with 
an "educated guess" for the nucleon propagator which can be improved perturbatively 
as higher powers of 5 are considered. Since one is looking for the propagator that leaves 
the energy stationary, one is not allowed to substitute the nucleon equation of motion 
in the expression of the energy. We could also interpolate kinetic energy terms, as in 
Ref. but we prefer to start with the simplest interpolation. Such generalizations 
will be examined in forthcoming publications. 

The fact that the divergent part of a physical quantity P{^) calculated with this 
interpolated Lagrangian will be 5 and /i dependent implies that the coefficients appear- 
ing in the counterterm Lagrangian Cct{^) will also depend on both parameters. The 
interpolated Cct{S) has the same field-operator structure as the original Cqt of Eq. (^, 
with 6- and yU-dependent coefficients. The explicit 6 and fi dependence is not important 
since this dependence will appear automatically in the process of fixing the value of the 
coefficients in the renormalization process. 

Next the strategy consists in calculating the interacting propagator, which is obtained 
by inverting Dysons's equation S^^{p) = ^ (p) — where the self-energy is 

calculated as a perturbation expansion in powers of 6, using for the "non-interacting" 
propagator, S^, the one corresponding to Cq. 5° can be split in the usual way as 
S^{k) = Sp{k) + S'l){k), where Sp{k) and Sl){k) are the Feynman and density dependent 
parts for quasi-particles of mass Vl and energy Eq = {1? + fi^)^/^. Inversion of Dyson's 
equation is standard 0. First, one notes that the in- medium is of the general 
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form 0]: S(A;) = S*(A;) —'y^T,^{k) +7-kS^(A;). Then, one defines the auxihary quantities: 



n*{k) = n + j:'{k) 

E*{k) = \k*^ + Q*\k) 



1/2 



k* = k[l + S'"(A;)], 

k*" = F + W{k) = \k^ + S°(A;), k* 



The propagator can be written as S{k) = Spik) + S^ik), with: 



Spik) 



}t* + n*{k) 



m- 



5{k^ - E{k))e{kF- |k|) , (9) 



}i* -VL*{k) + ie ' ' E*{k) 

where E{k) is the single-particle energy, which satisfies E{k) = [E*{k) — S°(A;)]^o^^^^j. 
Note that we have assumed that the nucleon propagator has simple poles with unit 
residue. Within the approximation scheme we will work in this paper, this assumption 
is satisfied, as can be seen below. 

The renormalization procedure in the 6 expansion follows closely the usual perturba- 
tive renormalization program, with small, however important, differences. Hence, in the 
following we shall repeat some standard textbook material which we feel essential for the 
appreciation of the differences. We start with the problem in vacuum {kp = 0), in which 
case S^{k) = Sp{k). We perform our calculations in 2ct; = 4 — 2e dimensions using dimen- 
sional regularization techniques To 0{6), the self-energy is simply ^^^^{k) = —6jj,, 
which arises from the bilinear term 6 fiipip. To 0(5^), one has in principle tadpole and 
exchange contributions: 



^(2)tad 



(k) = 



m. 



Tr 



9v 



^,Ti\YS%{q)\, (10) 



and 



T.f-\k)=^5'glj ^^Sf {q)/^.{k - q) - ^5^gl j ^,l,Sf\q)^.{k - q)Y ■ (H) 



d\ 



(2vr) 



The vector meson tadpole contribution vanishes after taking the trace and performing 
the integral, whereas scalar- meson tadpole contains a finite and a divergent part. The 
renormalization of this is performed using the term from Cct{^) in Eq. (|^). Its 

contribution to the self-energy is S^^ = (jf^ As(0). In order to produce a stable 

vacuum, (3i{5) is chosen so as to eliminate the divergent and finite parts. Any finite 
piece must be canceled, since such a term would contribute to the energy of the vacuum, 
which in turn could be lowered without bound. 
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The exchange parts of the self-energy are given by: 



^(2)exch 



'121 



where 



9l 



(47r)2 
(47r)' 
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(47r)2 
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where Ca and C;, are irrelevant constants and 2 = 1,2 are the integrals: 



{da)i In 



A;2(a - l)a + (1 - a)Q? + «m^) 



(13) 



(14) 



where ((ia)o = cxda and ((ia;)i = (ia , with a being the Feynman parameter. 

The renormalization of the exchange self-energy is carried out with the last two terms 
of Eq, (H), which contribute to the self-energy as S^'^(A;) = -M^{5) + C,i{^){^ - M), and 
the renormalized vacuum contribution to the self-energy, Tip{k), up to 0{6^) is given by 
j:R(k) = -6fi + sg)<=^^\A;) + SCT(A;). 

Generally, counterterms are composed of a divergent part which completely elim- 
inates the poles and of an arbitrary finite part which is fixed according to a chosen 
renormalization scheme |jll[. Here we shall consider the parameters gm,^m and M to be 
"the" renormalized physical couplings and masses. This choice amounts to the on-mass 
shell renormalization scheme in which the counterterms remove both divergent and finite 
contributions from loop corrections to measurable amplitudes. Within this renormal- 
ization scheme the finite parts of both counterterms are fixed by the renormalization 
conditions: 

dSp\k) 







and 



1. 



(15) 



Application of these conditions leads to the following expression for the renormalized 
self-energy in vacuum: 

Y^f{k) = A{k^)}i- B{k^), (16) 
where the functions A{k'^) and B{k'^) are free of divergencies and given by: 



A{k'' 



-6' 



aiP) - 2Mc{M^) 



(17) 



with a{P) = a{k^) - a{M^), b{P) = b{P) - b{M^), and c{M^) = Ma'{M^) - b'{M^), 
where the prime denotes a derivative with respect to the argument. It is easy to check 
that Eqs. (0) are indeed satisfied by Eqs. (|T6|-|1^. 

In principle one can proceed to higher orders in 6 and include, for example, the 
important vertex corrections. The renormalization at higher orders introduces no extra 
complications as compared to the usual perturbative renormalization. Instead of going to 
higher orders, we consider next the problem of the energy density, where new divergencies 
arise and extra renormalization is required. In order to simplify the discussion, we shall 
consider here only direct terms and then compare the results with the ones obtained 
within the RHA. 

When considering direct terms only, up to 0{6'^) one obtains the following relations 
for the auxihary quantities defined in Eq. (|): 



k 



0* 



ml Jo {2TT)^En{qy 



^ 2 / ^3' k* = k (19) 

m^j Jo [ZTry 

From these, one constructs the nucleon propagator which is then used in the expression 
for the energy density. We start with the baryon contribution to the energy density: 
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27r2 
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fi^F d^k P + MVt* 



B .fB 
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(20) 



lo (27r)3 E*{k) ^"^^ ' 
Note that the divergent and finite parts are different from the usual ones because we 
have not used the nucleon field equation of motion. Carrying out the VEV subtraction, 
the remaining divergencies are eliminated by the counterterm J2n=2iPn/n\)(j)"' in Eq. (^. 
Thus, the finite contribution to the baryonic part of the energy density is given by (after 
taking S = 1): 

8^ 



kF d^k k2 + M^]* 



(27r)3 E*{k) 
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where 



47r2 [^""{m) 
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+ M^{n* - M) - ^M^{n* - Mf - yM((]* - Mf - ^{n* - Mf 



[22) 



Proceeding in the same way for the scalar meson contribution, we obtain: 



2gl \ mlJ {2'nfE*{k) 



-Ac 



(23) 



where 



A 



9l 1 



n* 

M 



- M\Q* - M) - -M(Q* - Mf -—{pr- Mf 
2 6 



(24) 



Finally, there are no divergencies with the vector meson contribution, it is given by: 



1^ 

2 ml 



2 r 9 n2 



(25) 



The determination of the unknown parameter yU follows from the PMS applied to the 
energy density: 

dS dS dn* dS , , 

^ ^ = 0, (26) 



d\x dfl* djji dQ* 
where we used the fact that dQ*/dfi ^ 0. Differentiation of £■ with respect to f2* leads 
to the following self-consistent equation for the effective nucleon mass i7*: 

^2 ^ d^k k2 dAc 



-^gl r d^k n* 



ml J (27r)3 E*{k) 



Ac 



4^^ 

m2 J (27r)3 E*%k) dVl* 



0. (27) 



:281 



To arrive at this result we made use of the following identity: 

dQ* ^ ' gl dQ* 

Clearly, Eq. (|27|) admits two solutions. The one that follows from the vanishing term in 
curly braces corresponds to the usual RHA self-consistent solution. It is not difficult to 
show that when this is substituted into the equation for the energy density, one obtains 
the RHA energy density: 

'-kl\\-' 



9l 



2m?, VStt^ 



VF 



(29) 



where lS.EyF is the energy density corresponding to vacuum fluctuations, given by: 



A£. 



VF 



n*^\n ( — j + M-'{n*-M) 



47r2 

7 25 ■ 

+ -M\Q* - Mf + —{Q* - M)"^ + —{Q* ~ M)" 
2 3 12 



(30) 
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With the parameter set [|, gl = 62.89, = 79.78, m, = 550 MeV, = 783 MeV and 
M = 939 MeV, nuclear matter is saturated at kp = 1.42 fm~^, with a binding energy of 
16 MeV per nucleon. 

The other solution to Eq. (l27|) is: 



44/^^. (31) 



dn* ml J (27r)3 E*^{k) 
Note that the solutions of this equation are independent of the parameters of the model; 
the effective nucleon mass Q* depends only on M and kp. The numerical solution of 
this equation gives an effective nucleon mass that is an increasing function of the nuclear 
density. Moreover, when this solution is substituted into the original expression for the 
energy density, and using the same parameters as the ones of the RHA, we find a new 
saturation point for nuclear matter at kp ~ 5.3 fm^^, which corresponds to a density of 
approximately 50 times the one of normal nuclear matter, with a binding energy of the 
order of 570 MeV per nucleon. Of course, if the parameters are changed, these numbers 
will change. However, we do not pursue the discussion on this further because exchange 
terms should be included for consistency (they are of the same order in 6 as the direct 
ones), although one should notice that in all known cases where exchange graphs are 
included their effect can be absorbed by a readjustment of parameters 0. 

When one considers the exchange diagrams in the energy density, the renormalization 
becomes more difficult, but not more complicated than the usual renormalization in the 



loop expansion []T2|. In particular, if instead of inverting Dyson's equation for obtaining 
the propagator from the self-energy one uses the same approach as Chin [|l^] for treating 
the exchange diagrams, the renormalization is indeed simplified. 

We have given the first step towards the implementation of the (5-expansion for study- 
ing vacuum effects in effective hadron field theories. We demonstrated how the renor- 
malization program can be implemented within the 6 expansion. We also calculated 
the renormalized energy density of nuclear matter at 0{6^) by neglecting exchange dia- 
grams, and found that the usual RHA approximation is readily reproduced. In addition, 
a new state of nuclear matter is found from the PMS applied to the renormalized en- 
ergy density. Concluding, we believe that the results are very encouraging for future 
investigations concerning the vacuum of hadronic models. 
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